Gravitational Wave Production At The End Of Inflation 
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We consider gravitational wave production due to parametric resonance at the end of inflation, or 
"preheating" . This leads to large inhomogeneities which source a stochastic background of gravita- 
tional waves at scales inside the comoving Hubble horizon at the end of inflation. We confirm that 
the present amplitude of these gravitational waves need not depend on the inflationary energy scale. 
We analyze an explicit model where the inflationary energy scale is ~ 10 9 GeV, yielding a signal 
close to the sensitivity of Advanced L1GO and BBO. This signal highlights the possibility of a new 
observational "window" into inflationary physics, and provides significant motivation for searches 
for stochastic backgrounds of gravitational waves in the Hz to GHz range, with an amplitude on the 
order of Q gw (k)h 2 ~ 10~ n . Finally, the strategy used in our numerical computations is applicable 
to the gravitational waves generated by many inhomogeneous processes in the early universe. 



A successful model of inflation must have a "grace- 
ful exit" that describes the transition from the accel- 
erated phase to a thermalizcd universe Q. A widely 
studied mechanism for achieving this is preheating (e.g. 

U M M M M a ! S, M EH El El El El EH). Af- 

tcr inflation, the inflaton (or a related field) oscillates 
about the bottom of its potential, driving the resonant 
amplification of specific momentum modes of some cou- 
pled field(s). This renders the universe inhomogeneous, 
and the resulting spatial gradients source gravitational 
waves. For GUT inflation, the present peak frequency is 
between 1 MHz and 1 GHz 0, El- It was conjectured 
that the characteristic frequency is inversely proportional 
to the inflationary scale, while the amplitude can be in- 
dependent of this scale (IBj . leading to a signal poten- 
tially detectable by future iterations of LIGO and BBO. 
We confirm this conjecture by numerically computing the 
gravitational wave spectrum in a toy model of preheating 
following low scale inflation. While simple inflationary 
models typically involve GUT scale physics, many stringy 
models have a much lower inflationary scale, so this signal 
may eventually lead to new constraints on these models. 
The tools developed for this analysis will allow us to ex- 
plore fully realistic preheating models in an expanding 
background. Furthermore, our computational strategy 
applies to any inhomogeneous phase in the universe, and 
may have applications beyond the present problem. 

Computational Strategy & Results: During para- 
metric resonance, momentum modes of a field x are 
pumped by an oscillating field <p. In simple models <p 
is the inflaton, but in hybrid models 4> is the direction 
orthogonal to the inflationary trajectory which induces 
the "waterfall" transition 2C], l2l| . In either case the la- 
grangian can be expressed as 

c = - l -(d4>) 2 - l -{d x ?-v{<p lX ). (i) 

We numerically simulate the nonlinear field evolution in 
a conformally rigid spacetime background. We can then 
compute the spatial parts of T^ v at any given time. The 



tensor contribution to the metric perturbation and Ein- 
stein equations read 

ds 2 = dt 2 - a 2 (t) (Sij + hij) dx l dxi ,(2) 
G„„ (t) + 5G^ (x,t) = SnG [% v (t) + ST^ (x,t)] (3) 

where the overbar denotes the homogeneous background 
values. The perturbation is Transverse-Traceless, so 
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The Fourier components of the obey 
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and keeping only the transverse-traceless source. 

We evolve the fields with LatticeEasy [22| , treating 
the spatial background as a rigid, expanding box. We 
assume the fields' initial inhomogeneity is derived from 
their quantum mechanical vacuum fluctuations. The 
obey ordinary differential equations, which we solve 
with a fourth order Runge-Kutta integrator. Unlike 
UEH, we use no flat-space results, and the gravitational 
wave amplitude can be followed throughout the simula- 
tions. Numerical noise is a significant problem, since the 
transverse-traceless source involves differences of terms 
involving derivatives, which are themselves computed via 
numerical differencing on the lattice - which can easily 
lead to a loss of numerical significance. In future work 
we plan to perform simulations at higher spatial reso- 
lution, with different algorithms for the field evolution, 
in order to produce highly accurate spectra for a vari- 
ety of models over a large frequency range. Moreover, 
we have made simplifying assumptions about the sub- 
sequent evolution of the universe, which can rescale the 
overall tensor spectrum. In models without resonance, 
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FIG. 1: We plot the spectrum of gravitational radiation produced during resonance with 11 = 10~ 18 (left) through to 10" 6 
(right) in units where mpi ~ W 19 GeV = 1, where each spectrum has a value of /i 10 3 times large than the one immediately to 
the left. The corresponding initial energy densities run from from (4.5 x 10 9 GeV) 4 to (4.5 x 10 15 GeV) 4 for our choice of 4>o- 
The plots are made on 128 3 grids, and the "feature" at high frequency is a numerical artifact. 



the computed tensor spectrum is many orders of mag- 
nitude smaller than seen in resonant models, providing 
a "null check" of our algorithm, and we have recovered 
the key results of [H,[l5[. Compared to these results, the 
spectra in [H, [l5| are noticeably noisy, which is an ar- 
tifact of the older algorithm. The amplitudes obtained 
with our new code are somewhat lower than the previous 
results, possibly due to better noise suppression. 
We calculate the power spectrum using (23j 
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Evolving the h%j and inserting (JT]) into equation (20) of 
[l5| (and correcting a typographical error) we find: 



vt gw {k)h z = n r h z [ — 
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Here a e is the scale factor at the end of the simulation, fi r 
is the present-day density of radiation, go and g± denote 
the number of thermal degrees of freedom at present and 
at matter-radiation equality, and h is the dimensionless 
Hubble constant at the present epoch. We thus assume 
that the universe is radiation dominated between the end 
of our simulation and matter-radiation equality, during 
which time £l gw (k) is constant, and matter dominated 
thereafter. In the future we plan to consider more real- 
istic transfer functions [3, H^]. Moreover, we implicitly 
assume that the universe promptly thermalizes after pre- 
heating, which requires a trilinear coupling we have not 
included here 
We work with 



v(^x) = \n 2 4> 2 + \g 2 4> 2 x 2 - 



(9) 



In standard quadratic chaotic inflation, /j, is fixed by the 
power spectrum, and is not a free parameter. Since we 



are interested in the end of inflation, equation @ need 
not describe the potential in the epoch when cosmological 
perturbations were generated, and /x is a free parameter. 
To see this, consider 
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During inflation, <j> is large and a = 0. When <f> = M/h, 
the field rapidly evolves towards the minimum at a = 
M/y/X. This is the "waterfall" phase transition. Further, 
we assume that a = (a) and that there is there is no a 
production via a a 2 \ 2 term, giving 
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This is a mass term for (f>. In practice, we need m <C M 
in order to get the correct perturbation spectrum, so the 
effective mass is simply y? = h 2 M 2 /X. Here, we will 
simply take /j to be a free parameter. The inflationary 
energy density is then approximately p 2 (f> 2 l , where <fio is 
the field value at the beginning of our simulation, which 
we assume to be roughly coincident with the end of in- 
flation. We take 0o ~ -2m p /, so the inflationary energy 
scale is set by y}l 2 . For hybrid inflation, M/h ~ .2m p i 
and for fi = hM/VX = lCr 18 m p ; - fOGeV we need 
h 4 /X «3x 10~ 16 , which demands a small but not minis- 
cule value of h, if we assume A ~ 0.1. 

We hold the resonance parameter q = g 2 m 2 l /n 2 fixed 
[i~5| . This ensures that the same modes (in units of post- 
inflationary energy scale) are in resonance as \x is varied. 
We choose q = 2 x 10 6 , which makes g unrealistically 
small as \i is decreased. However, our purpose is to give 
an "existence proof" that (p)reheating can generate a 
substantial gravitational wave background at low infla- 
tionary scales. This assumption ensures that the struc- 
ture of resonance does not change with the inflationary 
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energy scale, allowing us to isolate the aspects of gravi- 
tational wave production which depend only on the over- 
all inflationary scale. Figure [T] shows the gravitational 
wave spectra for different inflationary energy scales. As 
predicted [15j, the peak frequency scales inversely with 
energy, whereas the amplitude is scale independent. 

The spectrum declines very rapidly at large k. This 
reflects the structure of the resonance bands, as modes 
with large k are never in resonance. Conversely, we see a 
relatively broad plateau, possibly attributable to modes 
moving through the resonance band as inflation contin- 
ues. The amplitude of modes which are super-horizon 
throughout resonance is expected to rise as k 3 or faster 
[ill ]. The universe grows significantly during resonance, 
so we cannot simultaneously resolve both resonant and 
super-horizon modes. In order to prevent contamination 
by numerical noise we must choose our box to ensure that 
we resolve a significant number of high-frequency modes 
which are never resonantly amplified. A very large scale 
simulation is thus needed to see the full shape of the peak. 

The source term is potentially non-zero whenever the 
fields are inhomogenous, but its amplitude rises rapidly 
during resonance, and then decreases as the universe be- 
gins to thermalize. Most of the gravitational waves are 
thus radiated during a comparatively narrow time in- 
terval. In the models treated here, the universe grows 
~ 50% larger as 90% of the power is generated. 

Discussion: Bubble-wall collisions after a first order 
phase transition at the end of inflation 
can also generate a present-day Q gw (k)h 2 ^ 
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spectrum has similar properties to the preheating back- 
ground: the frequency scales inversely with the inflation- 
ary scale, while the amplitude need not depend on the in- 
flationary scale. This process yields a maximum £l gw (k) 
similar to that found here for parametric resonance, al- 
though the detailed spectra may differ considerably. This 
is to be expected since, since gradient terms source hij 
and the amplitude of these is bounded by their contribu- 
tion to the total energy density. If the gradient energy 
saturates this bound, VL gw (k)h? ~ 10~ n or 10~ 10 is a 
gen eric result, in the absence of non-perturbative effects 
[33L [34I ]. In our simulations the gradient energy peaks 
at around 25% of the total energy density. Since there 
are necessarily kinetic and potential contributions to the 
energy, this fraction can not be significantly increased. 

Since our algorithm works for an arbitrary extended 
source of gravitational radiation, we could apply it to 
other mechanisms that generate stochastic back gro unds. 
These include networks of cosmic superstrings [301 ] . the 
TeV scale phase transition present in Randall-Sundrum 
models [31|, a first order electroweak phase transition 
[3^ ] . cosmological bubble collisions (28jj, or large-scale 
turbulence in the presence of magnetic fields [H, 0, HH . 
Moreover, |37[ presents an alternative (and physically 



distinct) mechanism for generating gravitational waves 
at the end of inflation, and we expect that it will be 



possible to analyze this mechanism with our algorithm. 
Finally, [361 ] discusses the generation of gravitation waves 
by density perturbations during the radiation dominated 
era, and it would be interesting to explore the overlap be- 
tween this second order calculation and our techniques. 

The detectability of high frequency background of 
stochastic gravitational radiation is an open question. 
We plot a schematic version of our results in Figure[2] and 
we see that the proposed space-based detectors BBO and 
LISA arc sensitive to much longer wavelengths than any 
signal likely to be generated during preheating from GUT 
scale inflation. It will be a stretch for Advanced LIGO 
to see the signals computed here. However, a further it- 
eration of LIGO is likely to put significant constraints 
on any signal that would be generated during resonance 
after inflation occurring at scales around 10 9 GeV. 

Detecting high frequency gravitational waves is partic- 
ularly challenging, since the required strain-sensitivity at 
fixed Q gw (k) scales as the cube of the frequency. Conse- 
quently, it is easier to detect a background arising from 
preheating after low-scale inflation. With few exceptions, 
stringy models of inflation have V 1 / 4 < 10 13 GeV, two or- 
ders of magnitude below GUT scales, and can be much 
lower [3^. Moreover, reheating after thermal inflation 
would lead to signal very close to the LIGO band [l7j |. 
In these models the primordial tensor spectrum is un- 
observable, but any preheating signal is easier to detect. 
Very low scale inflation (at or near the TeV scale e.g. 
[H Gil) could produce a signal visible to a BBO-style 
mission. Finally, future detector technologies might be 
sensitive to high-frequency gravitational radiation. 

Our key result is an existence proof that parametric 
resonance following inflation can lead to a significant 
gravitational wave background, independently of the in- 
flationary energy scale. Many questions remain. Firstly, 
we have only considered models which are well-described 
by {9]), and our set-up is a toy model, requiring a very 
small coupling parameter at low inflationary scales. This 
ensures the structure of resonance is independent of the 
energy scale, which allows us to isolate the impact on 
the inflationary energy scale on the height and location 
of the peak. However, this small coupling is not a prereq- 
uisite for preheating or resonance. Preheating is gener- 
ically associated with enhanced inhomogeneities which 
will source gravitational radiation, so we can hope that 
similar results also apply to more realistic resonance sce- 
narios, but these need to be explored in detail. Secondly, 
we have not yet considered overlapping backgrounds from 
other astrophysical processes. Thirdly, we need to incor- 
porate the finer details of the tensor mode transfer func- 
tion and the post-resonance thermalization of the uni- 
verse. 
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FIG. 2: We sketch the gravitational wave spectra obtained for the lowest and highest energy models computed here, relative 
to that of the Advanced LIGO goal, and the proposed LISA and BBO experiments. We see that inflationary models with lower 
energy scale may lead to a signal which is visible at LIGO scales if the sensitivity of LIGO is further improved, and with BBO. 
The tensor background generated by quantum fluctuations during GUT scale inflation is shown by the solid horizontal line. 
The dashed lines denote the inferred fe 3 tails. The spectra generated by the inflationary scenarios considered in 0, [l5[ roughly 
overlap with the 10 15 GeV spectrum depicted above. 
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